Magnetic excitations in two-leg spin 1/2 ladders: experiment and theory 
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Magnetic excitations in two- leg S=l/2 ladders are studied both experimentally and theoretically. 
Experimentally, we report on the reflectivity, the transmission and the optical conductivity o(tJ) of 
undoped La I Cai4_ :z ,Cu2404i for x=4, 5, and 5.2. Using two different theoretical approaches (Jordan- 
Wigner fermions and perturbation theory), we calculate the dispersion of the elementary triplets, 
the optical conductivity and the momentum-resolved spectral density of two-triplet excitations for 
0.2 < J\\/Jj_ < 1.2. We discuss phonon-assisted two-triplet absorption, the existence of two-triplet 
bound states, the two-triplet continuum, and the size of the exchange parameters. 



I. Introduction 

Low-dimensional quantum antiferromagnets offer a di- 
verse view on quantum fluctuations at work. Partic- 
ular interest has focused on two-leg spin 1/2 ladders, 
which show a spin liquid ground state with a spin gap 
to the lowest excited state, and superconductivity un- 
der pressure upon hole doping In undoped two- 
leg ladders the existence of two- "magnon" (two-triplet) 
bound states was predicted theoretically by a number of 
groups H-[|. Recently, we have reported on the obser- 
vation of a two-triplet bound state with St o t=0 in the 
optical conductivity spectrum of (La,Ca)i4Cu2404i ||, 
a compound with stacked layers of CU2O3 ladders and of 
CuC>2 chains. Concerning the peak positions and the line 
shape of the bound states, good agreement was achieved 
between experiment and two different theoretical ap- 
proaches, namely, Jordan- Wigner fermions and pertur- 
bation expansion using unitary transformations [^). Us- 
ing the second approach, some of us have discussed the 
momentum-resolved spectral densities and the relation- 
ship between fractional (spinons) and integer excitations 
(triplets) in two- leg spin 1/2 ladders pj. In these pro- 
ceedings, we discuss recent experimental and theoretical 
progress that we have achieved. 

II. Experiment 

Single crystals of La x Cai4_ x Cu2404i were grown by 
the travelling solvent floating zone method These 
so called telephone number compounds have attracted 
particular interest due to the possibility of hole dop- 
ing. Here, we are interested in the magnetic properties 
of nominally undoped samples, i.e. Cu 2+ , which corre- 
sponds to x=6. Single phase crystals could only be syn- 
thesized for x < 5.2 11 1 . The samples studied here with 
£=5.2, 5, and 4 on average contain rt=0.8/24, 1/24 and 



2/24 holes per Cu, respectively. However, x-ray absorp- 
tion data show that for these low doping levels the holes 
are located within the chains [Q. Thus, we consider the 
ladders to be undoped . 

In order to determine the optical conductivity <j(u>) 
we have measured both transmission and reflection data 
between 500 and 12 000 cm -1 on a Fourier spectrometer. 
The transmittance T{oj) and the reflectivity R(to) were 
calibrated against a reference aperture and a gold mirror, 
respectively. The optical conductivity o~(cu) = nKu/2Tr 
was determined by inverting [|l3|jl4| 

R(u) = [(n - l) 2 + K 2 ] I [{n + lf + K 2 ] , (1) 
T(u;)=[(l-R)H] I [l-(i?$) 2 ] , (2) 



$(w) = cxp(— 2und/c) = exp(— ad) , (3) 

where n denotes the index of refraction, k the extinction 
coefficient, a the absorption coefficient, c the velocity 
of light and d the thickness of the transmission sample 
[-R(w) denotes the single bounce reflectivity and hence 
needs to be measured on a thick ("semi-infinite"), opaque 
sample]. Equation 2 is obtained for a sample with par- 
allel surfaces by adding up the intensities of all multiply 
reflected beams incoherently, i.e., by neglecting interfer- 
ence effects. Experimentally, this condition is realized 
either if the sample surfaces are not perfectly parallel 
or by smoothing out the Fabry-Perot interference fringes 
by means of Fourier filtering. In case of weak absorp- 
tion k <C n, the reflectivity is entirely determined by n 
and not suitable to derive k by using a Kramers-Kronig 
transformation. At the same time, k can be determined 
very accurately from the transmission. Since T{u>) de- 
pends exponentially on k ■ d, the appropriate choice of d 
is essential. Furthermore, the thickness d determines the 
effective "resolution" given by the interference fringes. 
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The period of the fringes is given by Av = (2nd) -1 , 
where Av is given in wave numbers if d is given in units 
of cm. Thus, a large value of d has the advantage of nar- 
row fringes with a small amplitude, but strongly restricts 
the maximum value of k that can be observed. 

The top panel of Fig. [I] shows the reflectivity mea- 
sured on a 0.8 mm thick sample (x=A) for polarization 
of the electrical field both parallel to the rungs and to the 
legs. At about 600-700 cm" 1 one can observe the oxygen 
bond stretching phonon mode. At higher frequencies, 
the reflectivity is featureless, which is characteristic for 
the weak absorption regime below the gap of an insu- 
lator. We plotted the data for both 4 K (thick lines) 
and 300 K (thin dashed lines), but in R(u>) they are al- 
most indistinguishable. The different absolute values of 
the two polarization directions reflect the difference in 
n, namely, n a « 2.3 and n c ~ 2.6. The weak absorp- 
tion features we are looking for can only be detected in 
a transmission experiment, which also reveals a strong 
temperature dependence. The middle and bottom pan- 
els of Fig. [l] show T(uj) measured on thin single crystals 
with x=A (d=60/im) and x—5.2 (d=28/im). The data 
in the middle panel were measured with a resolution of 
4 cm" 1 and show only small interference fringes, most 
probably because the sample surfaces were not perfectly 
parallel. The sample with x=5.2 (bottom panel) showed 
strong interference fringes with a period of « 70 cm , 
which have been removed by Fourier filtering. 

In order to obtain a qualitative estimate of the absorp- 
tion features it is common to approximate Eq. 2 by 
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FIG. 2. Bottom panel: a\[ui) of La a! Cai4_ ;!; Cu24 04i (x=4, 
5, and 5.2) at 4 K, calculated from transmission and reflectiv- 
ity after Eq. 2. Top and middle panel: Plotting the logarithm 
of the transmission divided by the sample thickness already 
gives a qualitative estimate of the absorption features without 
knowledge of the reflectivity. 



T(uj) « (1 -i?) 2 $ 



(4) 



In this case one obtains the absorption coefficient a(to) 
as 



a(u>) « - ln(T) jd + 2 m(l - R) /d 



(5) 



Neglecting the almost constant second term, we plot 
-ln(T(w))/dfor x=4, 5, and 5.2 at 4 K in the top (E \\ c) 
and middle panels (E \\ a) of Fig. E| For comparison, the 
optical conductivity er(w) as derived by inverting Eq. 2 
without further assumptions is displayed in the bottom 
panel of Fig. ||. In a(u>) the curves for x—5 and x=5.2 
almost fall on top of each other. Note that the erroneous 
discrepancy between these two data sets in — \n(T)/d is 
not due to a difference in R (u) , but due to the approxima- 
tion used in Eq. 4. A precise and reliable determination 
of o~(u>) in the case of small absorption thus requires the 
measurement of both the transmission T and the reflec- 
tivity R and the use of Eq. 2. 

III. Magnetic contribution to <t(u) 

The spectra of (t{uS) can be divided into three different 
regimes. The absorption below w 1300 cm" 1 can be at- 
tributed to phonons and multi-phonon bands. The strong 
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FIG. 3. Optical conductivity of La5.2Ca8.sCu24 04i 

at 4 K (thick lines) on a linear (top panel) and on a logarith- 
mic scale (bottom panel). The thin lines show Gaussian fits of 
the electronic background. For E || a the fit is unambiguous. 
For E [| c two different background fits are shown in order to 
illustrate the uncertainty (see Fig. W). 



increase at high frequencies is due to an electronic back- 
ground, which most probably has to be identified with 
the onset of charge-transfer excitations. We focus on the 
features in the intermediate frequency range. In order to 
analyze them we used Gaussian fits |L5| to subtract the 
electronic background (thin lines in Fig. 3|). For E \\ a, the 
background can be determined unambiguously by fitting 
the measured data for w>7000 cm -1 . After subtraction 
of the background the a a (u) curves are nearly indepen- 
dent of x, which corroborates the assumption that the 
ladders are undoped. For E \\ c, the higher absorption 
complicates the determination of the background con- 
siderably. We were able to measure the transmission 
up to « 6200 cm -1 (in case of the x=5.2 sample with 
d=28 /im). Although this is a significant improvement 
compared to our earlier work which was restricted to 
u> < 5000 cm -1 (a;=5), it is not yet sufficient for an un- 
ambiguous determination of the background. Therefore 
we show two different possibilities for Gaussian fits of the 
c-axis background in Fig. |[ 

The estimates of the magnetic contribution to cr(u> 
which we have obtained this way are plotted in Fig 
(thin and thick solid lines). In Ref. |J we interpreted 
these features in terms of phonon-assisted two-triplet ab- 
sorption g(J, which was used to describe a\ of the un- 
doped 2D cuprates (e.g. YBa 2 Cu 3 6 @) and of the ID 
S=l/2 chain S^CuOs [Gl|. Because of spin conserva- 
tion two triplets are excited. The simultaneous excita- 
tion of a phonon provides the symmetry breaking neces- 
sary to bypass the selection rule and at the same time 
takes care of momentum conservation [[) 16 l^] . This 
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FIG. 4. Comparison of the magnetic contribution to the 
optical conductivity of La5.2Ca8.sCu2404i at 4 K (solid lines) 
and of the calculated results obtained with Jordan- Wigner 
fermions {J\\/ J± = 1.0; Jm =1100 cm -1 ; dashed) and with op- 
timized perturbation theory (Jy/Ji = 1.15; Jm = 1080 cm -1 ; 
dotted). In the top panel, the two different experimental 
curves (thin and thick solid lines) correspond to the two al- 
ternative background fits of Fig. 0. 

process is the lowest order magnetic response possible 
in infrared absorption. The two peaks between 2000 and 
3000 cm -1 (see Fig. |j) can be identified with the ID van 
Hove singularities in the density of states of the strongly 
dispersing two-triplet bound state with S to t—0 |J (see 
below). The absorption at higher energies is attributed 
to the two-triplet continuum. The improved estimate 
of the background as compared to Ref. II allows us to 
detect two peaks within the continuum at about 3800 
and 4700 cm -1 for E \\ c. The precise determination 
of the spectral weight of these high energy features re- 
quires additional work. Note that also the a-axis spec- 
trum displays a shoulder at high energies (sa 5000 cm -1 ), 
although it is quite weak. 

IV. Theory 

Antiferromagnetic S —1/2 two-leg Heisenberg ladders 
are represented by the Hamiltonian 



jj| H — { J|| (Si^Si^+i + S2,iS2,,:+i) + J_i_Si,iS2,i 



}, (6) 



where Jj_ and Jm denote the rung and leg couplings, re- 
spectively. For Jii=0 one can excite local rung singlets 
to rung triplets which become dispersive on finite Jii (see 
Fig. ||). For Ji=0 the S=l chain excitations decay into 
asymptotically free 5=1/2 spinons. An intuitive picture 
of the "magnons" (elementary triplets) for J±, J|| ^ 
can be obtained from both limits: the elementary exci- 
tations are either dressed triplet excitations or pairs of 
bound spinons with a finite gap A as long as Jj_ > p0) . 
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FIG. 5. Elementary triplet dispersion for < J\\/J± < 1.2. 
For each theory, the thick lines denote the coupling ratio 
that offers the best description of the optical conductivity 
of (La,0a)i4Cu24O4i (see Fig. 0). 



Here, we present two theoretical approaches to describe 
the excitations of the ladder. One approach makes use 
of the Jordan- Wigner transformation to rewrite the spins 
as fermions with a long-ranged phase factor. Expanding 
the phase factor yields new interaction terms between 
the fermions. The resulting interacting fermion problem 
is treated diagrammatically. A similar treatment works 
very well for a ID chain. The other approach uses an ex- 
trapolated perturbation in J\\/J±_ and separates contri- 
butions of different triplet number by continuous unitary 
transformations |8|JlCt|. Both methods are controlled in 
the sense that they become exact on Jm / J± — » 0. In fact, 
comparing the results of these two methods for the dis- 
persion of the elementary triplets (see Fig. |), the agree- 
ment is excellent for J\\/J± < 0.6. 

In order to determine the optical conductivity one 
needs to calculate the momentum-resolved two-particle 
spectral densities with S t ot = The evolution of 

the spectral densities from J\\/J± = 0.2 to J\\/J± = 1 
is plotted for both theories in Fig. ^. The optical con- 
ductivity as given in Fig. |^ is obtained by integrating 
these fc-resolved curves with a weight factor oj ■ sin (fc/2) 
For comparison with experiment one has to add the 
phonon energy as a constant shift of the energy scale 
(before multiplication by to). 

For J\\/J± = 0.2 the elementary triplet still shows only 
little dispersion (see Fig. |J). Therefore the two-particle 
continuum is rather narrow; the spectral weight is piled 
up at the bottom of the continuum for small momenta 
k, and for large k a bound state is formed below the 
continuum (see Fig. [?]). For small ratios of J\\/J± this 
bound state reaches its maximum energy at the Brillouin 
zone boundary and dominates in a(to) as a single sharp 
peak with only a small continuum contribution at higher 
energies H (sec left panels of Fig. |). With increasing 
J\\/Jj_ the bound state acquires a strong dispersion, and 
for J\\/J± > 0.5 it shows a maximum at k « tt/2 and 



Jordan- Wigner - 

- 3 




2 12 3 

co ( J ± ) co ( J ± ) 

FIG. 6. Comparison of the optical conductivity for E || c 
(leg) calculated with Jordan- Wigner fermions (top pan- 
els) and with optimized perturbation theory (bottom) for 
< J\\/ J-L < 1.15. The left panels focus on the bound state, 
whereas the right panels emphasize the two-triplet continuum 
on an enlarged j/-scale. For clarity, the curves were shifted 
with respect to each other by 3 (0.5) (f2cm) _1 in the left 
(right) panels. For each theory, the thick lines denote the 
coupling ratio that offers the best description of the bound 
states in ct(lu) of (La,Ca)i4Cu2404i (see Fig. |ij). For compar- 
ison with experiment, the phonon frequency still has to be 
added (here: 600 cm' 1 ). 

a minimum at the Brillouin zone boundary (see Fig. [?]). 
Both give rise to van Hove singularities in the density of 
states which cause peaks in <t(w). Therefore, the domi- 
nant peak observed for J\\/J±_ = 0.2 splits into two with 
increasing J|| (see left panels of Fig. ||). Comparison of 
these spectra with the experimental data offers a very ac- 
curate tool to determine the exchange coupling constants, 
since the frequencies of the two peaks depend strongly 
on the coupling ratio J\\/J± ||- For La 5 . 2 Ca8.8Cu2404i 
we obtain J\\/J± = 1.0 and J11 = 1100 cm -1 from the 
Jordan- Wigner fermions, whereas perturbation theory 
yields J\\/J± = 1.15 and Jy = 1080 cm" 1 . This 15% 
discrepancy reflects the differences between the two the- 
ories that were already present in the dispersion of the 
elementary triplet (see Fig. ||). Both theories reproduce 
the experimental line shape rather well for both polariza- 
tion directions. Note that the lower peak is suppressed 
for E || a due to symmetry 

A ratio of J\\/J± ~ 1 seems to be in conflict with 
several former results of other techniques, proposing 
J\\/J±. ^1-5 (see discussion in Q). Such a rather small 



4 



value of J± is suggested by the small spin gap observed in 
e.g. neutron scattering |2C ] . On the basis of our results we 
can exclude J\\/J± > 1.2 [E|. Recently, it was pointed out 
that the neutron data (i.e., the small spin gap) are also 
consistent with an isotropic exchange J\\/J±_ ~ 1 — 1.1 
and J|| s» 900 cm -1 , if a ring exchange of J cyc s» 0.15 
J|| is taken into account |20|. Adding a finite ring ex- 
change term reduces the gap at k = ir and washes out the 
dip in the triplet dispersion at small k 20 2l]]. Further- 
more, the cyclic exchange weakens the attractive inter- 
action between two rung triplets. We expect that these 
two changes, reduction of the dip and of the attractive 
interaction, render the theoretical predictions closer to 
the experimental findings of optical spectroscopy |^l|,22|. 
Since a « 15% ring exchange term is estimated [P0]2l| 
we expect a similar change in Jy/Jx. 

Finally, we address the evolution of the continuum with 
increasing Ju/Jx (see Fig. 0). At small k, the spec- 
tral density is broadened strongly by the increase of the 
continuum width. Therefore, the k—0 Raman response 
shows a sharp peak for small J\\/J± and a broad band for 
J|i/J_L«1 (see Ref. [Q for more details). For large mo- 
menta one can observe the opposite, the features within 
the continuum become stronger and more pronounced 
with increasing J\\/J±. The mid-band line of square root 
singularities which runs from k « 7r/2 to k — n denotes 
the upper edge of the precursor of the spinon contin- 
uum that is well known from the chain limit JlO| . For 
J|i/Jj_=l the perturbation result (top left panel of Fig. 

shows additional pronounced features within the high 
energy part of the continuum JTo[ ] . The appearance of 
these features is related to the existence of the dip in the 
dispersion of the elementary triplet at small k (Fig. ||). 
Precursors of these features are present in perturbation 
theory for J||/Jj_=0.8, where the dip in the dispersion 
is only small. The Jordan- Wigner fermions do not show 
these features. 

In the weighted superpositions (Fig. |^) the perturba- 
tive results for J\\/J± > 0.8 display a second peak above 
lu « 3 J± in the high energy continua. The experimental 
data (Fig. |j) also display two peaks in the continuum, 
independent of the precise background correction. Yet 
the experimental peaks lie at higher energies and they 
are further apart from each other. On inclusion of ring 
exchange the calculated peaks will be shifted to higher 
energies since the attractive interaction is weakened. So 
we tentatively identify the two experimental with the two 
theoretical peaks even though their shape and mutual 
distance are not in perfect agreement. Further investiga- 
tions including ring exchange are indispensable. 

Summarizing, we reported in detail on one of the first 
observations of bound states in gapped spin systems. By 
phonon-assisted two-triplet absorption we were able to 
detect a bound state in the S = channel resulting 
from two elementary triplets in the undoped spin lad- 
der La x Cai4_ x Cu2404x. The bound state exists for mo- 



menta k 2l 0.37T. Thus, it was decisive that phonons are 
involved in the optical excitation process in order to pro- 
vide the necessary momentum. Two theoretical methods 
were employed to trace the evolution of the spectral den- 
sities as function of the coupling ratio Jn/J±. In both 
approaches the overall position and weight of the bound 
state agree well with the measurements. In experiment 
and in the perturbative approach two peaks in the con- 
tinua are found. Some quantitative differences indicate 
that an extension of a simple spin ladder Hamiltonian is 
necessary. We argued that a 10-15% ring exchange term 
is the appropriate extension. 

This project is supported by the DFG (FR 754/2-1, 
SP 1073 and SFB 484), by the BMBF (13N6918/1) and 
by the DAAD within the scope of PROCOPE. 
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optimized perturbation Jordan- Wigner fermions 
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FIG. 7. Momentum dependence of the two-triplet spectral densities with Stot = for E || c. Calculations were performed with 
optimized perturbation (left) and with Jordan- Wigner fermions (right) for < J\\/J± < 1. The fc-resolved and u-integrated 
weights for the two theories agree within 25% (note the different scalings of the plots) . Gray curves represent the bound state, 
which was divided by 16 and artificially broadened by Jx/100 in the left panels. Black curves indicate the continuum. Inset of 
bottom left panel: enlarged view of the continuum, the bound state is not shown here. 
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